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ABSTRACT 

We present a detailed strong lensing, weak lensing and X-ray analysis of Abell 2744 
(z = 0.308), one of the most actively merging galaxy clusters known. It appears to 
have unleashed 'dark', 'ghost', 'bullet' and 'stripped' substructures, each ~ 1O 14 M . 
The phenomenology is complex and will present a challenge for numerical simulations 
to reproduce. With new, multiband HST imaging, we identify 34 strongly-lensed im- 
ages of 11 galaxies around the massive Southern 'core'. Combining this with weak 
lensing data from HST, VLT and Subaru, we produce the most detailed mass map 
of this cluster to date. We also perform an independent analysis of archival Chan- 
dra X-ray imaging. Our analyses support a recent claim that the Southern core and 
Northwestern substructure are post-merger and exhibit morphology similar to the 
Bullet Cluster viewed from an angle. From the separation between X-ray emitting gas 
and lensing mass in the Southern core, we derive a new and independent constraint 
on the self-interaction cross section of dark matter particles cr/m < 3 ± 1 cm 2 /g. In 
the Northwestern substructure, the gas, dark matter, and galaxy components have 
become separated by much larger distances. Most curiously, the 'ghost' clump (pri- 
marily gas) leads the 'dark' clump (primarily dark matter) by more than 150 kpc. We 
propose an enhanced 'ram-pressure slingshot' scenario which may have yielded this 
reversal of components with such a large separation, but needs further confirmation 
by follow-up observations and numerical simulations. A secondary merger involves a 
second 'bullet' clump in the North and an extremely 'stripped' clump to the West. 
The latter appears to exhibit the largest separation between dark matter and X-ray 
emitting baryons detected to date in our sky. 

Key words: gravitational lensing: strong — gravitational lensing: weak — galaxies: 
clusters: individual: Abell 2744 — dark matter — large-scale structure of the Universe 
— X-rays: individual: Abell 2744. 
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The standard ACDM cosmological model suggests a bottom- 
up sequence of structure formation, in which a series of 
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merging events culminates in massive clusters of galaxies, 
the latest structures to form in the observable Universe 
(Bond et al. 1991; Lacey & Cole 1993). The number of clus- 
ters as a function of their mass (the steep, high end of the 
mass function) depends sensitively upon cosmological pa- 
rameters (e.g. Vikhlinin et al. 2009) and has become an im- 
portant observational test of cosmology. Most measurements 
of cluster masses rely upon the calibration of more easily ob- 
servable proxies, such as X-ray luminosity, temperature or 
galaxy richness. However, clusters form through multiple, 
dynamic accretions, so are likely to be turbulent places, and 
the turmoil affects those observable proxies. It is therefore 
vital to quantitatively understand the merging process, for 
example by mapping the distribution of dark matter, stars 
and baryonic gas in systems at many different stages of the 
merger. 

Merging clusters of galaxies have become useful labo- 
ratories in which to study the nature and interaction prop- 
erties of dark matter. The best-studied example is the Bul- 
let Cluster 1ES 0657-558 (z = 0.296) (Tucker et al. 1998; 
Markevitch et al. 2002). Combined X-ray and gravitational 
lensing analyses show a clear separation between the centres 
of X-ray emission and the peaks in surface mass density - 
indicating a fundamental difference between baryonic gas, 
which feels the pressure of the collision, and dark matter, 
which is nearly collisionless (Clowe et al. 2004, 2006; Bradac 
et al. 2006). The discovery and interpretation of the Bul- 
let Cluster has inspired a lively debate about whether such 
a system could exist in different cosmological models (see 
e.g. Hayashi & White 2006). Improvements are continuing 
in numerical simulations (Milosavljevic et al. 2007; Springel 
& Farrar 2007; Mastropietro & Burkert 2008; Lee & Ko- 
matsu 2010). Observations have also broadened, with discov- 
eries of a possible line-of-sight merger CL0024+1654 (Czoske 
et al. 2002; Hoekstra 2007; Jee et al. 2007; Zitrin et al. 
2009; Umetsu et al. 2010; Zu Hone et al. 2009, a), and other 
systems including the Baby Bullet (MACS J0025.4-1222) 
(Bradac et al. 2008), the Cosmic Train Wreck (Abell 520) 
(Mahdavi et al. 2007; Okabe & Umetsu 2008), Abell 2146 
(Russell et al. 2010), Abell 521 (Giacintucci et al. 2008; Ok- 
abe et al. 2010a) and Abell 3667 (Finoguenov et al. 2010). 
All these systems place potentially tight constraints on the 
interaction between baryons and dark matter, and are ex- 
emplary probes for our understanding of structure formation 
within gravitationally bound systems. 

In the archival data of 38 merging clusters, Shan et al. 
(2010) found the largest offset between X-ray and lensing 
signals to occur in the massive (Lx = 3.1xl0 45 erg/s in the 
2-10 keV range, Allen (1998)) cluster Abell 2744 (also known 
as AC118, or RXCJ0014.3-3022) at a redshift of z = 0.308 
(Couch & Newell 1984). The - 250 kpc offset is larger than 
that in the Bullet Cluster - although, as we shall discuss 
later, this value does not describe the separation of the 
main mass-clump from its stripped gas component. Nev- 
ertheless, Abell 2744 is undergoing a particularly interest- 
ing merger. The complex interplay between multiple dark 
matter and baryonic components appears to have unleashed 
'ghost', 'dark', 'stripped' and 'bullet' clusters. 

The first hint that Abell 2744 is in the middle of a ma- 
jor merging event arose from observation of a powerful and 
extended radio halo (P(1.4GHz) > 1.6 x 10 36 Watt, Gio- 
vannini et al. (1999); Govoni et al. (2001, a)). This indicated 



the presence of relativistic electrons accelerated through 
high Mach shocks or turbulence (e.g. Sarazin 2004). The 
picture was clarified by X-ray studies (Kempner & David 
2004; Zhang et al. 2004) that revealed substructure near the 
cluster core, plus an additional luminous structure towards 
the Northwest. Kinematic observations of cluster member 
galaxies (Girardi & Mezzetti 2001) suggested a bimodal dis- 
tribution in redshift space, but were not at first consid- 
ered significant. Recent kinematic studies focussing solely 
on Abell 2744 definitely show a bimodal velocity dispersion 
in the cluster centre, together with a third group of cluster 
members near the Northwestern X-ray peak (Boschin et al. 
2006; Braglia et al. 2009). Although there is no evidence for 
non-thermal X-ray emission (Million & Allen 2009), the frac- 
tion of blue star-forming galaxies (Braglia et al. 2007, and 
references therein) also seems to be enhanced. A default ex- 
planation emerged for the centre of Abell 2744, featuring a 
major merger in the North-South direction with a small in- 
clination towards the line-of-sight and a ^ 3:1 mass ratio of 
the merging entities (Kempner & David 2004; Boschin et al. 
2006). More controversial is the role of the Northwestern 
structure. Kempner & David (2004) detected a cold front 
on its SW edge and a possible shock front towards the clus- 
ter core, so concluded that it is falling towards the main 
mass. More recent analysis of Chandra data (Owers et al. 
2011) failed to confirm the presence of the shock front, only 
a cold front towards the Northern edge, and the authors 
proposed that the structure is moving instead towards the 
North/Northeast, after being deected from the main cluster 
in an off-centre core passage. What seems sure is that we 
are at least observing a complicated merger between three 
separate bodies (Braglia et al. 2007). 

Until now, Abell 2744 has been better constrained from 
X-ray and kinematic studies than by gravitational lensing. 
So far, Smail et al. (1997) detected a weak- lensing signal and 
strong-lensing features, followed by Allen (1998) who found 
a large discrepancy in the mass estimates for Abell 2744 
from X-ray and strong-lensing reconstructions. Given the 
cluster's dynamical state, this finding is perhaps no longer 
surprising since the merging would induce non-thermal sup- 
port and elongation along the line-of-sight, which increases 
the systematics from both methods. The most recent weak- 
lensing analysis (Cypriano et al. 2004) showed indications of 
substructure in the reconstructed surface-mass density, but 
did not reach the resolution required for more quantitative 
statements. Here, we present the results of an HST imag- 
ing survey, aimed at clarifying the evolutionary stage of this 
complex system. 

This article is organised as follows. In Sec. 2 we present 
our comprehensive lensing analysis, mainly based on newly 
acquired data taken with the Advanced Camera for Surveys 
(ACS) on the Hubble Space Telescope (HST). In Sec. 3 we 
describe a complementary X-ray analysis based on Chandra 
data. We discuss cosmological implications and an interpre- 
tation of the merging scenario in Sec. 4, and we conclude 
in Sec. 5. Throughout this paper we assume a cosmologi- 
cal model with Qm — 0.3, Qa — 0.7 and h — 0.7. At the 
cluster's redshift z — 0.308, one arcsecond corresponds to 
4.536 kpc. 
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2 LENSING ANALYSIS 

Our recently acquired, multiband HST/ACS imaging en- 
ables us to significantly improve upon previous mass mod- 
els of Abell 2744. Using the parametric method described in 
Sec. 2.2, we have identified strong gravitational lensing of 11 
background galaxies producing 34 multiple images around 
the Southern core, with an Einstein radius of ~ 30" (see 
below) . These enable us to tightly constrain the position and 
shape of the core mass distribution. No such multiple image 
systems are revealed around the N or NW clumps, immedi- 
ately indicating that their masses are lower. Our HST im- 
ages also yield ~ 62 galaxies/arcmin 2 for weak lensing anal- 
ysis (after charge transfer inefficiency (CTI) corrections are 
performed, as described below). This enables detailed mass 
modelling throughout our HST field of view. To probe the 
cluster merger on even larger scales, we incorporate ground- 
based weak lensing measurements from VLT and Subaru. 

We simultaneously fit all of these strong- and weak- 
lensing observations using our well-tested mass reconstruc- 
tion algorithm (Merten et al. 2009; Meneghetti et al. 2010) 
(which is similar to that used by Bradac et al. (2006) and 
Bradac et al. (2008) to map the Bullet Cluster and the Baby 
Bullet). Importantly, we make no assumptions about mass 
tracing light in the combined analysis. It should be noted 
though, that we rely on a parametric method, assuming that 
light traces mass, to identify multiple image systems. We 
account for possible uncertainties in this identification, by 
resampling the set of less confident identifications in the er- 
rors analysis based on a bootstrapping technique (see Sec. 
2.4). Our analysis reveals four individual clumps of mass 
1O 14 M0 within a 250 kpc radius. Previous weak lensing 
analysis of VLT images alone had resolved but a single broad 
mass clump (Cypriano et al. 2004). Below we describe our 
datasets, analyses, and results in more detail. The central 
cluster field and a preview on the matter and gas distribu- 
tion is presented in Fig. 1, which also shows the contours 
encircling different areas of specific probability, that a mass 
peak is located at this specific location of the reconstructed 
field. 



2.1 The HST/ACS dataset 

The HST data consist of two pointings in Cycle 17 
(data taken between Oct. 27-30 2009, Proposal ID: 11689, 
P.L: R. Dupke) with ~ 50% overlap between the point- 
ings. The images were taken with the ACS/WFC camera 
using three different filters, F435W (16.2 ksec 1 ), F606W 
(13.3 ksec) and F814W (13.2 ksec). 

The HST/ACS camera had been in orbit for eight years 
when the imaging was acquired. During this time above 
the protection of the Earth's atmosphere, its CCD detec- 
tors had been irreparably damaged by a bombardment of 
high energy particles. During CCD readout, photoelectrons 
are transported to the readout amplifier through a silicon 
lattice. Damage to this lattice creates charge traps that de- 
lay some electrons and spuriously trail the image - in a way 
that alters the shapes of galaxies more than the gravitational 
lensing signal that we are trying to measure. To undo this 
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trailing and correct the _raw images pixel-by-pixel, we used 
the detector readout model of Massey et al. (2010), updated 
for device performance post Servicing Mission 4 by Massey 
(2010). The corrected data was then reduced via the stan- 
dard CALACS pipeline (Pavlovsky 2006), and stacked using 
multidrizzle (Koekemoer et al. 2002). 

2.2 Strong lensing 

We concentrate our strong-lensing analysis on the reduced 
ACS images only. Several strong-lensing features are imme- 
diately identifiable by eye on the combined three-band image 
(Fig. 2). To find additional multiple images across the field 
of view, we apply the well-tested approach of Zitrin et al. 
(2009) (see also Zitrin et al. 2011b) to lens modelling, which 
has previously uncovered large numbers of multiply-lensed 
galaxies in ACS images of Abell 1689, C10024, 12 high-z 
MACS clusters, MS 1358 and Abell 383 (respectively, Broad- 
hurst et al. 2005; Zitrin et al. 2009, 2011, 2010, 2011a,c) . 

In the Zitrin et al. (2009) method, the large-scale distri- 
bution of cluster mass is approximated by assigning a power- 
law mass profile to each galaxy, the sum of which is then 
smoothed. The degree of smoothing (S) and the index of 
the power-law (q) are the most important free parameters 
determining the mass profile. A worthwhile improvement in 
fitting the location of the lensed images is generally found 
by expanding to first order the gravitational potential of this 
smooth component, equivalent to a coherent shear describ- 
ing the overall matter ellipticity, where the direction of the 
shear and its amplitude are free parameters. This allows for 
some flexibility in the relation between the distribution of 
dark matter and the distribution of galaxies, which cannot 
be expected to trace each other in detail. The total deflec- 
tion field C£t(0), consists of the galaxy component, a ga i(0), 
scaled by a factor K ga i , the cluster dark matter component 
(Xdm(0), scaled by (1-K ga i), and the external shear compo- 
nent a ex (0): 

a T (0) = K gal a gal 0) + (1 - K gal )a DM {0) + a ex (0), (1) 

where the deflection field at position m due to the external 
shear, a ex (6 m ) = (a eXjX , a ex , y ), is given by: 

OLex,x(0m) = \j\ cos(20 7 )Ax m + \j\ sin(20 7 ) Ay m , (2) 

(9 m ) = |7l sin(20 7 )Ax m - | 7 | cos(20 7 )Ay m , (3) 

where (Ax m , Ay m ) is the displacement vector of the posi- 
tion 6m with respect to a fiducial reference position, which 
we take as the lower- left pixel position (1,1), and 7 is the 
position angle of the spin- 2 external gravitational shear mea- 
sured anti-clockwise from the x-axis. The normalisation of 
the model and the relative scaling of the smooth dark matter 
component versus the galaxy contribution brings the total 
number of free parameters in the model to 6. This approach 
to strong lensing is sufficient to accurately predict the lo- 
cations and internal structure of multiple images, since in 
practice the number of multiple images readily exceeds the 
number of free parameters, which become fully constrained. 

Two of the 6 free parameters, namely the galaxy power 
law index q and the smoothing degree S, can be initially set 
to reasonable values so that only 4 of the free parameters 
need to be fitted at first. This sets a very reliable starting- 



4 J. Merten et al 






Figure 1. The field of Abell 2744, with different interesting features. The false-colour background is provided by HST/ACS (the two 
pointings can be identified by the higher, bluely background noise level), VLT and Subaru images on a field size of 240 ;/ x 240 ;/ (~ 1.1 
Mpc on a side). Overlaid in cyan are the surface-mass density contours most concentrated in the 'core' area and in magenta the more 
evenly-distributed X-ray luminosity contours. The peak positions of the core, N, NW, and W clumps are indicated by the green likelihood 
contours, derived from the bootstrap samples. Contours are 86%, 61%, and 37% of the peak likelihood for each clump. (Assuming a 
Gaussian probability distribution, these would correspond to 0.3, 1, and 2-sigma confidence contours.) Note that the NW clump most 
likely peaks in the area indicated by NW1 (in 95% of the bootstrap realisations), but peaks also at NW2 in 54% of the bootstrap 
realisations. The small red circles show the position of the local overdensities in the gas distribution, associated with each individual dark 
matter clump. The white rulers show the separation between dark matter peaks and the bright clusters galaxies and local gas peaks. 



point using obvious systems. The mass distribution is there- 
fore well-constrained and uncovers many multiple-images 
that can be iteratively incorporated into the model, by using 
their redshift estimation, from photometry, spectroscopy or 
model prediction and location in the image-plane. At each 
stage of the iteration, we use the model to lens the most ob- 
vious lensed galaxies back to the source plane by subtracting 
the derived deflection field, then relens the source plane to 
predict the detailed appearance and location of additional 
counter images, which may then be identified in the data by 
morphology, internal structure and colour. We stress that 
multiple images found this way must be accurately repro- 
duced by our model and are not simply eyeball "candidates" 
requiring redshift verification. The best fit is assessed by the 
minimum RMS uncertainty in the image plane 



^ M5 images = X^fo ~ x ^ + ^ ~ y ^ I ^images' ( 4 ) 

i 

where and y i are the locations given by the model, xi 
and yi are the real image locations, and the sum is over 
the total number of images -7Vj ma g es - The best-fit solution 
is unique in this context, and the model uncertainty is de- 
termined by the locations of predicted images in the image 
plane. Importantly, this image-plane minimisation does not 
suffer from the well known bias involved with source plane 
minimisation, where solutions are biased by minimal scatter 
towards shallow mass profiles with correspondingly higher 
magnification. 

In Abell 2744 we have uncovered a total number of 34 
multiple images belonging to 11 background sources. We 
label the different systems in Fig. 2, which also shows the 
critical curve of the cluster derived from the strong lensing 
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Table 1. The multiple-image system of Abell 2744. 



Figure 2. A zoom into the innermost core region of the 
HST/ACS images. Shown as a continuous white line is the critical 
curve of the cluster as it is derived from the strong-lensing model. 
It assumes a source redshift zs = 2.0. Also shown are the approx- 
imate positions of the identified multiple-image systems as they 
are listed in Tab. 1 with a varying two-colour scheme to avoid 
confusion of close- by images. These systems were found by the 
method of Zitrin et al. (2009) and are not simple optical identi- 
fications. The visible field size is ~ 100" x 100 /; , translating to 
~ 450 kpc on a side. 



model. The model predicts an Einstein radius of = 30" ± 
3". With such a large number of clear multiple images we 
can constrain the inner mass distribution very well, and we 
shall incorporate these constraints into our joint strong- and 
weak-lensing analysis described in Sec. 2.4. 



2.3 Weak lensing 

Several areas outside the cluster core are of special interest, 
due to the complicated structure of this merging cluster. 
We measure the shapes of background galaxies to derive a 
weak-lensing signal and extend our mass reconstruction over 
a much larger field. Since our HST data cover only a limited 
field of view, we also include VLT and Subaru imaging in 
our weak lensing analysis. Additional HST pointings in the 
future are crucial for a full interpretation of this merger, as 
we shall discuss in the course of this analysis. 



2.3.1 HST/ACS 

To select lensed background galaxies, we obtained photo- 
metric redshifts based on our three niters using Bayesian 
photometric redshifts (BPZ) (Bemtez 2000; Coe et al. 2006). 
Cluster ellipticals at z ~ 0.3 occupy a unique region in our 
colour-colour space, enabling us to obtain better than ex- 
pected results. 

Of 118 galaxies with published spectroscopic redshifts 
(Owers et al. 2011 and references therein; all z < 0.7) within 
our field of view, 99 yield confident photo-zs, accurate to 





Image-ID 


X 


77 


z 


v\ so1 


irce) . (image^ (additional knot)) 


\ ) 


\ ) 






1.1 


-35.10 


-13.55 


2.0 ±0.3 




1.2 


—30.15 


—23.95 






1.3 


0.10 


-35.40 






1.11 


-33.60 


-16.50 






1.21 


-31.50 


-21.55 






1.31 


1.60 


-35.85 






2.1 


9.30 


-11.50 


2.0 ± 0.3 




2.2 


-34.25 


12.35 






2.3 


2.80 


1.05 






(2.4) 


-0.50 


-7.10 






2.11 


11.55 


-7.65 






2.21 


-32.55 


13.90 






2.31 


5.55 


3.10 






(2.41) 


-0.05 


-4.15 






3.1 


-10.10 


22.65 


4.0 ± 0.3 




3.2 


-7.40 


22.95 






(3.3) 


27.15 


2.20 






4.1 


-18.25 


-9.00 


3.5 ± 0.3 




4.2 


-29.25 


-5.30 






4.3 


18.05 


-31.60 






5.1* 


8.85 


29.00 


4.0 ± 0.5 




5.2* 


3.85 


31.60 






5.3* 


19.65 


19.30 






6.1* 


-38.15 


-5.95 


3.0 ± 0.5 




6.2* 


-24.25 


-28.30 






6.3* 


-0.60 


-33.15 






7.1* 


-37.35 


-7.85 


3.7 ± 0.5 




7.2* 


-27.85 


-26.10 






7.3* 


5.10 


-34.85 






8.1 


— 10.70 


20.90 


4.0 ± 0.2 




8.2 


-8.00 


21.40 






(8.3) 


16.10 


5.90 






9.1* 


-6.65 


-18.45 


3.0 ±0.5 




9.2* 


-2.80 


-21.90 






(9.3) 


-43.20 


10.80 






10.1* 


-6.65 


-20.65 


3.0 ±0.5 




10.2* 


-3.55 


-22.75 






10.3* 


-44.90 


11.00 






(11.1) 


-16.00 


-13.30 


3.0 ±0.5 




(11.2) 


-34.20 


-4.65 






(11.3) 


10.70 


-31.55 





If the image ID is shown in brackets, the image is not confi- 
dently reproduced by the lensing model described in Sec. 2.2. 
The systems marked with an asterisk, were nicely reproduced by 
the parametric lensing model but we allow them to be resam- 
pled by our nonparametric reconstruction technique to allow for 
uncertainties in the number of identified strong lensing features. 
All other images define a 'confident catalogue' of multiple-image 
systems and are used in our subsequent analysis. The x-and y- 
coordinates are relative to the BCG position («j2000 = 3.58611°, 



°J2000 



-30.40024°) in arcseconds. Redshifts of each system 



and their respective error are derived from the model predictions. 
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Az rsj 0.06(1 + z) RMS with no significant outliers. Back- 
ground galaxies were selected as those with confident photo- 
z > 0.5. This cut successfully excludes all 5 foreground 
and 86 cluster galaxies in our spec-z sample with confident 
photo-z (though one cluster galaxy was assigned a less con- 
fident photo-z 0.8). However, given our limited photomet- 
ric coverage, we do not rule out some contamination of our 
background sample with foreground and/or cluster galaxies. 

We measure the weak gravitational lensing signal in the 
F814W HST exposures, which are the deepest and contain 
the most galaxies at high redshift. We measure their shapes, 
and correct them for convolution by the Point- Spread Func- 
tion (PSF), using the 'RRG' (Rhodes et al. 2000) pipeline 
developed for the HST COSMOS survey (Leauthaud et al. 
2007; Massey et al. 2007). The RRG method is particularly 
optimised for use on high resolution, space-based data. Since 
HST expands and contracts as it warms in the sun or passes 
through the shadow of the Earth, even this telescope does 
not have a constant PSF. However, 97% of the variation in 
its PSF can be accounted for by variation in its focal length 
(Jee et al. 2007, the separation between the primary and 
secondary mirrors). We therefore measure its focal length 
by matching the shapes of the ~ 12 bright stars in each 
pointing to models created by raytracing through the opti- 
cal design (Krist 2003) . This achieves a repeatable precision 
of 1/xra in the determination of the focal length, and we con- 
struct a PSF model from all those stars observed during the 
600-orbit COSMOS survey at a similar focus (Rhodes et al. 
2007). We finally use this PSF model to correct the shapes 
of galaxies, and to obtain estimates of the amount by which 
their light has been sheared. The result is a catalogue of 1205 
galaxies with shear estimates, corresponding to a density of 
~ 62 galaxies/arcmin 2 . 



2.3.2 VLT/FORS1 

The complementary VLT data for our weak lensing analysis 
is identical to that included by Cypriano et al. (2004) in a 
study of 24 X-ray Abell clusters. The total field of view is 6.8 
on a side, centred on the brightest cluster galaxy (BCG) and 
significantly exceeding the coverage of our HST imaging. V, 
R and /-band imaging was obtained with the FORS1 camera 
between April and July 2001, with exposure times of 330 s in 
each filter. The data were reduced with standard IRAF rou- 
tines and, to maximise the depth, we perform weak- lensing 
shape analysis on the combined VRI image. Seeing condi- 
tions were excellent, with a FWHM of stars in the combined 
VRI image of 0759. 

We measure galaxy shapes and perform PSF correc- 
tion using the IM2SHAPE method (Bridle et al. 2002). This 
involves a two-step process to first map the PSF variation 
across the observed field using stars, then to model each de- 
tected galaxy, perform PSF correction and recover its ellip- 
ticity. To remove foreground contamination and unreliable 
shape measurements from our catalogue, we apply magni- 
tude cuts plus additional rejection criteria (see Cypriano 
et al. 2004). Since that work, we have improved the effi- 
ciency of foreground galaxy removal and now keep a higher 
density of background galaxies in our shear catalogue. The 
result is a catalogue of 912 galaxies with shear estimates, or 
~ 20 galaxies/arcmin 2 . 



2.3.3 Subaru/ SuprimeC am 

To extend the total field of view even further, especially 
in the Northern areas of the cluster field, we obtained 1.68 
ksec i'-band imaging data with Subaru/ SuprimeC am dur- 
ing Semester S08B. The data were reduced following Okabe 
& Umetsu (2008) and Okabe et al. (2010,a). Astrometric 
calibration was conducted by fitting the final stacked image 
with the 2MASS data point source catalogue; residual as- 
trometric errors were less than the CCD pixel size. Due to 
poor weather conditions, the seeing size is as large as l'/28. 

We measure galaxy shapes and perform PSF correction 
using the IMCAT package (provided by Kaiser et al. 1995) in 
the same pipeline as Okabe et al. (20 10, a) with some modi- 
fications following Erben et al. (2001). Background galax- 
ies were selected in the range of 22 < % < 26 ABmag 
and f£ + cr r * ~ 3.4 < ru < 6.0 pix, where ru is the 
half-light radius, and f£ and a r * are the median and stan- 
dard error of stellar half-light radii m* , corresponding to 
the half median width of the circularised PSF. The den- 
sity of background galaxies in our final shear catalogue is 
^15 galaxies/arcmin 2 . This is 30-50% of typical values from 
images obtained during normal weather conditions Okabe 
et al. (2010,a). 



2.4 Combined lensing reconstruction 

In order to combine the weak- and strong-lensing constraints 
in a consistent way, we use the joint lensing reconstruction 
algorithm described in Merten et al. (2009) (see also Bradac 
et al. 2005, 2009, for a similar approach). This method 
has been extensively tested in Meneghetti et al. (2010) and 
proved its capability to faithfully recover the cluster mass 
distribution over a broad range of scales. 

Our joint mass reconstruction is nonpar ametric, in the 
sense that it neither makes any a priori assumptions about 
the cluster's underlying mass distribution nor does it need 
to trace any light-emitting component in the observed field. 
However, as described in Sec. 2.2, we use a parametric 
method to identify the multiple image systems in the field. 
We reconstruct the cluster's lensing potential (its gravita- 
tional potential projected onto the plane of the sky) ip by 
combining measurements of the position of the critical line 
and the reduced shear. To do this, we divide the observed 
field into an adaptive mesh, which discretises all observed 
and reconstructed quantities. A statistical approach is cho- 
sen to combine our various measurements, by defining a 
multi-component x 2_ f unc tion that depends on the underly- 
ing lensing potential and a regularisation term R(ip) to pre- 
vent the reconstruction overfitting noise (see Merten et al. 
2009; Bradac et al. 2005) 



x 2 W = xwW + xsW + #W. 



(5) 



The weak-lensing term is defined by the expectation 
value of the complex reduced shear in each mesh position 
(e), which is obtained by averaging the measured ellipticities 
of all background galaxies within that grid cell 



2 / Z(*bM ^ r-i ( i \ 



- Z(z) K (rl>) 



where Z(z) is a cosmological weight factor as defined e.g. in 
Bartelmann & Schneider (2001), and 7 = is the 



The galaxy cluster merger Abell 2744 7 



shear of the lens, with the two components expressed in 
complex notation, k = dd*?p is the convergence, where 
the complex differential operator in the plane is defined as 
d := ( t^- + i with 61 and 62 being the two angular coor- 
dinates in the sky. The indices i, j indicate the discretisation 
of the input data and the lens properties, where we have to 
take into account the full x 2 -function because the averaging 
process of background galaxies might result in an overlap 
of neighbouring mesh points, expressed by the covariance 
matrix dj. 



The strong-lensing term is defined as 



|2\ 2 



2 (det^)) 2 fc _ ((1 - Z(s)ttM) 2 - \Z(zhW\' Jk (7 , 

where the index k labels all pixels in the reconstruction 
mesh, which are supposed to be part of the critical curve 
within the uncertainties cr s , given by the pixel size of the 
grid. At these points, the Jacobian determinant detA(ip) of 
the lens mapping must vanish. 

We iterate towards a best-fitting lens potential by min- 
imising the x 2_ function at each mesh position 

d X 2 W 1 



dipi 



with I G [0, 



(8) 



In practice, we achieve this by translating this operation into 
a linear system of equations and invoking a two-level itera- 
tion scheme (see Merten et al. 2009, and references therein). 

In this analysis, we use all strongly lensed multiple- 
image systems from the confident sample (Tab. 1) together 
with their derived redshifts, and combine weak lensing shear 
catalogues from all three telescopes (see Sec. 2.3) as de- 
scribed in Sec. 2.4.1. The combination of strong lensing data 
and the density of background galaxies allows for a recon- 
struction on a mesh of 72 x 72 pixels in the central region 
(corresponding to a pixel-scale of 8.4"/pix) and 36x36 pixels 
in the outskirts of the field (corresponding to a pixel-scale of 
16.7"/pix). Error estimates were produced by bootstrapping 
the redshift uncertainties of the strong-lensing constraints, 
and resampling those multiple image systems in Tab. 1, 
which are marked with an asterisk. The result are 500 dif- 
ferent bootstrap realisations of the refined cluster core. 150 
bootstrap realisations of the cluster outskirts were produced 
by bootstrapping the combined ellipticity catalogues. The 
number of bootstrap realisations is mainly constrained by 
runtime considerations. All error estimates are calculated 
from the scatter within the full bootstrap sample. If not 
stated differently, the given value reflects 68% confidence 
level. 



2.4-1 Combining the catalogues 

In order to combine the three different catalogues of ellip- 
ticity measurements we used the following strategy. Clearly, 
the Subaru catalogue delivered the largest field-of-view but 
it was derived from only single-band imaging and under bad 
seeing conditions. Therefore, we decided to limit its field 
size to a box of 600" x 600" around the centre of the cluster 
and to cut out the central ~ 400" x 400" part, which was 
sufficiently covered by HST/ACS and VLT exposures with 
better data quality. As a result, Subaru data only covers the 
outermost 200"on each side of the field. 
One might argue that the combination of HST/ACS and 



VLT ellipticities in the innermost centre of the cluster is 
problematic, but indeed the galaxy density is three times 
higher in the HST/ACS field and clearly dominating the 
adaptive-averaging process of the nonparametric reconstruc- 
tion algorithm that we used for the further analysis (com- 
pare Sec. 2.4 and Merten et al. (2009)). 
Another issue is the possible double-counting of galaxies, so 
we cross-correlated both catalogues with a correlation ra- 
dius of 2.5"and found 160 double-count candidates over the 
full HST/ACS field. Given the pixel size of the final recon- 
struction, this translates to an insignificant average double- 
count probability of less than one galaxy per pixel. Fur- 
thermore, the weighting scheme of the adaptive- averaging 
process was implemented such, that the highest weight of 
the VLT ellipticities was identical to the smallest weight of 
the HST/ACS ellipticities. Finally, the errors in the ellip- 
ticity measurement are treated in the joint reconstruction 
method in a purely statistical way by deriving the variance 
of the weight- averaged sample of ellipticities in each recon- 
struction pixel. 

To derive a physical surface mass density from the scaled 
lensing convergence one needs to know at least the mean 
redshift of the background galaxy population that was used 
to produce the ellipticity catalogues. Problems with the dif- 
ferent depths of the fields and therefore with the final mass 
analysis should not be a crucial issue for a relatively low- 
redshift cluster like Abell 2744. However, the redshift for 
determining the surface-mass density from the reconstructed 
convergence in the overlap area has been calculated as the 
galaxy-density weighted average of both, the VLT and the 
HST/ACS populations. The redshifts of the strong lensing 
features (compare Tab. 1) were included for the determina- 
tion of the core mass and the respective redshifts of the VLT 
and Subaru source distribution was used in the outskirts of 
the field. 

In order to test the effect of dilution in the outskirts of 
the field, we increased the ellipticity values for all Subaru 
background galaxies by 15% and repeated the reconstruc- 
tion. This test is necessary due to the single-band Subaru 
imaging. As it turns out, the difference in the reconstructed 
convergence is marginal since the ellipticity values are al- 
ready low with large scatter in this area of the field. How- 
ever, the effect was included in the determination of the error 
budget for the reconstructed total mass and mass profile. 

2.5 Reconstruction results 

We obtain a map of the lensing convergence across the field 
(proportional to the projected mass) by applying the Lapla- 
cian operator to the lensing potential on the adaptively re- 
fined mesh (see Fig. 3). We find a total mass (by assuming 
h = 0.7 here and further on) within a radius of 1.3 Mpc 
around the Core of M(r < 1.3Mpc) = 1.8 ± 0.4 x 1O 15 M , 
which is in good agreement with kinematically derived 
masses (Boschin et al. 2006). A mass determination within a 
field of (1300 x 750) kpc (to compare easily to the work on the 
Bullet Cluster of Bradac et al. (2006)) centred on the Core 
density peak yields M(1.73Mpc 2 ) = 7.4 ± 1.0 x 10 14 M o , 
rendering Abell 2744 comparable in mass or slightly less 
massive than the Bullet Cluster (Bradac et al. 2006). The 
overall radial convergence and mass profile can be found in 
Fig. 4. 
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Figure 3. The convergence map of the cluster field for a source 
redshift of zs = oo and a field size of 600" x 600", translating 
to ~ 2.7 Mpc on a side. The black contours start at ko = 0.14 
with a linear spacing of Ak = 0.047. The four white circles with 
labels indicate identified sub-clumps and the radius within which 
their mass is calculated. The radius of all four circles is 55.4" w 
250 kpc. 

Most interestingly, our gravitational lensing analysis re- 
solves four distinct sub- structures, indicated in Fig. 3 by the 
white circles. We label these substructures as Core, North- 
western (NW, later on dubbed as 'dark'), Western (W, later 
on dubbed as 'stripped') and Northern (N) structure. Please 
note that there is some indication of a possible double peak 
in the NW structure, as can be seen in Fig. 1 and also in the 
clearly visible extension of the structure towards the West 
in Fig. 3. The Core, NW and W clumps are clear detections 
in the surface-mass density distribution with 11 a, 4.9a and 
3.8cr significance over the background level, respectively. 2 
Somewhat fainter with 2.3cr significance is the N structure, 
but it clearly coincides with a prominent X-ray substructure 
found by Owers et al. (2011) and is therefore included in the 
further analysis (compare Sec. 3). The positions of the mass 
peaks and their local projected masses within 250 kpc are 
listed in Tab. 2 and shown in more detail in Fig. 1. The new 
HST/ACS images thus allow a striking improvement in our 
map of the mass distribution and reveal the distribution of 
dark matter sub- structure in great detail for the first time. 
We will refer to the individual mass clumps resolved here in 
our discussion of the X-ray analysis below. 



3 X-RAY ANALYSIS 

The most difficult part of interpreting merging clusters is de- 
termining the geometric configuration of the collision, such 

2 The background level was estimated in the following way: From 
the total convergence, shown in Fig. 3, a radial area of 75" around 
the Core, 70" around NW and 55" around N and W was cut out. 
From the remaining field the mean and variance in the conver- 
gence level was calculated. 



Table 2. Structures identified within our lensing reconstruction. 



Name x y M(r < 250kpc) 

(") (") (10 14 M ) 



Core 


-13^ 


_4+ 6 

4 -13 


2.24 ±0.55 


NW 


n±\l 


84+.1 5 


1.15 ±0.23 


W 


177±l 3 7 


-30i- 


1.11 ±0.28 


N 




170tl° 8 


0.86 ±0.22 



The x-and y-coordinates are provided in arcseconds, relative to 
the BCG position (a J20 00 = 3.58611°, <5 J2000 = -30.40024°). 
The 68% confidence limits on peak positions are derived from 500 
bootstrap realisations for the Core peak, and from the pixel size 
of the coarse weak lensing mesh in the reconstruction outskirts for 
the other mass peaks. Masses assume h = 0.7 and their 68% con- 
fidence limits are derived from bootstrap realisations as described 
in the text. 




radial distance from main clump centre [kpc] 



Figure 4. Radial mass profiles of the cluster, derived from our 
lensing analysis. Shown in blue (referring to the left y-axis) is 
the convergence profile for a source redshift z s = oo. The large 
uncertainty for small radii arises from uncertainty in the exact 
position of the cluster centre. The cumulative mass as a function 
of radial distance from the cluster centre is shown in red (referring 
to the right y-axis). The black vertical line indicates the distance 
of the outermost multiple images from the cluster centre. 

as its impact velocity, impact parameter and angle with re- 
spect to the plane of the sky (Markevitch et al. 2002). For 
this purpose, X-ray data becomes a crucial addition to lens- 
ing measurements. The location of any shock fronts are re- 
vealed in the temperature of the intracluster medium (ICM). 
Velocities can be inferred from the density and temperature 
of intracluster gas (if the merger axis is near the plane of 
the sky), or through direct Doppler measurements (if the 
merger axis is near the line of sight). 

3.1 Reduction of the Chandra data 

We reanalysed all existing Chandra data of Abell 2744 
(listed in Owers et al. (2011)), using CIA0 4.3 with the cali- 
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bration database CALDB 4.4.2. We cleaned the data using 
the standard procedure 3 and kept events with grades 0, 2, 
3, 4 and 6. We removed the ACIS particle background as 
prescribed for 'VFAINT' mode, and applied gain map cor- 
rection, together with Pulse Hight Amplitude (PHA) and 
pixel randomisation. Point sources were identified and re- 
moved, and the sky background to be subtracted from spec- 
tral fits was generated from Blank- Sky observations using 
the acis_bkgrnd_lookup script. We fit spectra using XSPEC 
V12.6.0q (Arnaud 1996) and we adopt Vapec thermal emis- 
sion model from atomic data in the companion Astrophys- 
ical Plasma Emission Database (Smith et al. 2001), allow- 
ing for the variation of several individual elemental abun- 
dances during the spectral fittings. Galactic photoelectric 
absorption is incorporated using the wabs model (Morrison 
& McCammon 1983). Spectral channels are grouped to have 
at least 20 counts/channel. Energy ranges are restricted to 
0.5-9.5 keV. Metal abundances are quoted relative to the 
solar photospheric values of Anders Ez Grevesse (1989), and 
the spectral fitting parameter errors are 1-a unless stated 
otherwise. 

There is a known reduction of quantum efficiency at en- 
ergies below 1 keV due to a build-up of molecular contami- 
nants on the optical blocking filter (or on the CCD chips) 4 . 
To prevent this from affecting our measurement of low en- 
ergy line abundances, hydrogen column density and overall 
gas temperature, we fix the column density to the cluster's 
nominal value of 1.6 xlO 20 cm -2 (Dickey & Lockman 1990). 
To be conservative, we present results from only those re- 
gions of the CCDs best-suited for velocity analysis (ACIS-I 
pointings 8477 and 8557); for an independent sanity check, 
we also repeat the analysis with ACIS-S pointing 2212. We 
exclude pointings 7712 and 7915 because the regions of in- 
terest cross multiple CCDs, and interchip gain fluctuations 
could introduce spurious systematic effects when measuring 
the redshifts from X-ray spectra. 

Since spectral models are (weakly) degenerate with clus- 
ter redshift, it is common for simultaneous fitting routines 
to get stuck in local % 2 minima before they reach the global 
best fit. We circumvent this via an iterative approach. We 
perform initial fits while varying the redshift values within 
reasonable ranges (via the command STEPPAR in XSPEC). 
We then fix the redshift, and refit full spectral models to in- 
fer gas temperature, metal abundances and normalisations. 
Subsequently, we use these best-fit values as inputs in a new 
fit with the redshift free to vary. This provides a more re- 
liable estimation of the error on the redshift measurement. 
We repeat the process until the best-fit redshift no longer 
changes between iterations. 

3.2 Previous X-ray observations & interpretations 

Abell 2744 shows an extremely disturbed X-ray morphology 
(compare Fig. 5). With 25 ksec of Chandra ACIS-S3 imag- 
ing, Kempner & David (2004) decided it is the aftermath 
of a N-S collision between similar mass proto-clusters, at 
Mach > 2.6. They also found tentative evidence for a cold 

3 http://cxc.harvard.edu/ciao/guides/acis_data.html 

4 fit tp : / /cxc . harvard . edu /cal /Links / Acis / 
acis / CaLprods / qeDeg/index. html 



front in the detached nearby NW 'interloper', which they 
guessed was falling into the main cluster. Owers et al. (2011) 
obtained a further 101 ksec of ACIS-I imaging. The deeper 
data revealed the 'Southern minor remnant core' (SMRC) 
to be colder (Tx ~ 7.5 keV) than its surroundings, with a 
high temperature region to the SE (Tx > 15 keV) that they 
interpreted as a shock front. They concluded the SMRC had 
been a low-mass bullet that has passed through the 'North- 
ern major remnant core' (NMRC), leaving central tidal de- 
bris (CTD). The pressure ratio of ~ 3:1 across the shock 
front corresponds to a sky-projected shock velocity of 2150 
km/s (for an average temperature Tx ~ 8.6 keV). Owers 
et al. (2011) reversed the Kempner & David (2004) model 
of the interloper, concluding that it came originally from the 
South, has already passed through the main body of the host 
cluster with a large impact parameter, and is now climbing 
out towards the N-NW. 

Owers et al. (2011) also obtained spectra of more than 
1200 galaxies with the Anglo- Australian Telescope multi- 
fibre AAOmega spectrograph, and confirmed the velocity 
bi-modality of the cluster galaxies. One component of galax- 
ies near the SMRC has a peculiar velocity of 2300 km/s; a 
separate component near the NMRC has a peculiar veloc- 
ity of about — 1600km/s and enhanced metal abundances 
(~ 0.5 Solar). Assuming that the Northern and Southern 
cores have the same velocities as the 'apparently' associated 
galaxy populations, they de-projected the velocity of this 
collision to Mach 3.31 or nearly 5000 km/s. 

An interesting prediction of this scenario is that the in- 
tracluster gas should show a strong radial velocity gradient 
of ~ 4000 km/s or A z ~ 0.014 from North to South. Since 
the intracluster medium is enriched with heavy elements, ra- 
dial velocity measurements could be carried out directly by 
measuring the Doppler shift of emission lines in the X-ray 
spectrum (Dupke & Bregman 2001a, b; Andersson & Made- 
jski 2004; Dupke et al. 2007) or through changes in line 
broadening due to turbulence (Inogamov & Sunyaev 2003; 
Sunyaev et al. 2003; Pawl et al. 2005). The former requires 
high photon counts within the spectral lines and excellent 
control of instrumental gain but could, in principle, be mea- 
sured with current X-ray spectrometers. The latter requires 
very high spectral resolution that should become available 
through the future ASTRO-H and IXO satellites (however, 
see Sanders et al. 2010, 2011). 

The existing data are not ideal for measurements of 
ICM velocity structure with high precision due to the varia- 
tion of gain expected from the analysis of multiple observa- 
tions with different pointings and in different epochs. How- 
ever, the expected radial velocity gradient of 4000km/s, pre- 
dicted from the optical work of Owers et al. (2011) is higher 
than the expected interchip and intrachip gain variations, 
so that it becomes possible to test, even if just for consis- 
tency, the proposed merger configuration, as we did in this 
work. Chandra has a good gain temporal stability (Grant 
2001) and we shall control for spatial variations by perform- 
ing resolved spectroscopy of multiple cluster regions using 
the same CCD location in different ACIS-I pointings. 

3.3 Velocity measurements 

We perform our velocity analysis twice: for the regions of 
interest defined by Owers et al. (2011), then for the sub- 
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structures prominent in our lensing mass maps. These are 
respectively indicated by green or blue circles in Fig. 5. 

We find temperature and metal abundance values for 
the NMRC, CTD and SMRC regions of 8.29±0.67 keV and 
0.55±0.18 Solar ( x „ =1.06 for 234 degrees of freedom); 
9.51±0.47 keV and 0.27±0.08 Solar ( X » =1.03 for 469 de- 
grees of freedom); and 7.98±0.76 keV and 1.17±0.51 Solar 
(Xu =1.01 for 143 degrees of freedom). These values are 
derived from simultaneous fits to pointings 8477, 8557 and 
2212. The values are consistent with those obtained through 
individual spectral fits for each instrument/observation (two 
independent ACIS-I and one ACIS-S). Our temperature and 
abundance measurements are thus consistent with the val- 
ues found by Owers et al. (2011). Given that the observa- 
tions were not tailored for velocity measurements (i.e. we 
cannot exclude temporal or inter-chip gain variations), we 
conservatively include in our velocity error bars secondary 
X 2 minima that are separated from the global minimum at 
less than 90% confidence in the velocity measurements. 

Our velocity measurements are consistent with the pres- 
ence of a gradient between the NMRC and SMRC, but in 
the opposite sense to that expected from Owers et al. (2011) 
interpretation. Using the ACIS-S, the CCD with the best 
spectral response, the velocity difference between NMRC 
and SMRC is found to be > 1500 km/s at 90% confidence. 

The observations (8477 and 8557) indicate a higher 
magnitude for the velocity gradient detected. However, the 
central values of the best-fit redshifts in each observation 
are found to be significantly discrepant, introducing larger 
uncertainties. The reason for this discrepancy is not clear, 
but it is possibly related to the observed focal plane tem- 
perature variation of about 0.9C between these observations 
(Catherine Grant Personal Communication). This explana- 
tion is supported by the marginally significant observed sys- 
tematic difference of ~ 10%-20% in the best-fit values of 
temperatures and abundances measured for the same re- 
gions in these different ACIS-I observations. This effect was 
also noted in Owers et al. (2011) (Owers, M. Personal Com- 
munication). Therefore, although our analysis is consistent 
with a high velocity gradient, the direction of the gradient 
does not corroborate the merger configuration suggested by 
Owers et al. (2011). 

Having the lensing mass reconstruction shown in the 
previous sections at hand we can study the line-of-sight 
gas velocity mapping for the regions of interest (Core and 
NW) more precisely. Using all three exposures, we mea- 
sure their gas temperatures and metal abundances to be 
Tx Core = 10.50±0.57 keV, A Co r e =0.28±0.11 Solar with a 
reduced chi-squared x^— 0.90 for 466 degrees of freedom 
and T X7VW = 10.22±0.54 keV, Aivw=0.40±0.09 Solar with 
Xi/=0.91 for 477 degrees of freedom. 

Despite their similar gas properties, these regions show 
a velocity gradient > 5200 km/s at near 90% confidence, 
even including a conservative (la) intra-chip gain variation 
of 1000 km/s (e.g. Dupke & Bregman 2006) per CCD in 
the error budget, in quadrature. A contour plot of the ve- 
locity difference is shown in Fig. 6. This result is consistent 
with the idea that the Southern Core mass is redshifted with 
respect to other structures. Nonetheless, these calculations 
should be taken with caution due to the uncertainties related 
to temporal gain variations and inter-chip gain variations 
between ACIS-I and ACIS-S. The same analysis using only 




Figure 5. X-ray image of Abell 2744, overlaid with our lensing 
mass reconstruction (magenta). The velocity gradient is maximal 
between three regions of interest named by Owers et al. (2011) 
and circled in green (Southern minor remnant core, Northern ma- 
jor remnant core, and Central tidal debris). A velocity gradient 
is also detected between the regions of interest from our lensing 
mass reconstruction, circled in blue (Core and NW) as well as the 
interloper. 
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Figure 6. Simultaneous fits to the redshifts of the Core re- 
gion and the NW substructure using two ACIS-I and one ACIS-S 
pointings. The two contours correspond to 68% (inner) and 90% 
(outer) confidence levels, and the diagonal line indicates equal 
redshifts to guide the eye. The Core mass clump appears signifi- 
cantly redshifted when compared to the NW mass clump. 



the ACIS-I observations shows the same velocity trend, but 
with a reduction in the absolute velocity difference (an up- 
per limit of <4400 km/s at the 90% confidence level). Deeper 
Chandra observations, specifically tailored for velocity mea- 
surements will be crucial to further reduce uncertainties and 
to better constrain the velocity structure of this cluster. 
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4 INTERPRETING THE MERGER 

Abell 2744 is undoubtedly undergoing a complicated merg- 
ing process on a large cluster scale. Progressively more 
detailed studies, culminating in our lensing reconstruction 
(Sec. 2) and X-ray analysis (Sec. 3), have only agreed that 
the merger is more complex than previously thought. For 
example, the likely explanation for the gas velocity gradient 
in the opposite direction to that expected by Owers et al. 
(2011) is that the NMRC is not the main cluster. Our lens- 
ing mass reconstruction shows that the deepest gravitational 
potential is by far the Southern 'Core' structure, which is 
roughly coincident with Owers et al.'s SMRC but slightly 
further offset from the Compact Core. We also find three 
separate mass concentrations to the North, Northwest and 
West. Overlaying the lensing mass reconstruction and X-ray 
emission in Fig. 1 reveals a complex picture of separations 
between the dark matter and baryonic components. To in- 
terpret the sequence of events that led up to this present 
state, we shall now tour the regions of interest, with more 
detailed discussions. 



4.1 Core, the massive clump 

According to our lensing analysis, the Core region (lower- 
left quadrant of Fig. 1) is by far the most massive structure 
within the merging system (c.f. Tab. 2). All the strong lens- 
ing features can be seen within this clump. We find no large 
separation between the distribution of mass and baryonic 
components. The mass peak is centred amongst the bright 
cluster member galaxies (within la errors, it is consistent 
with the position of the BCG) and only 22 ± 12 arcseconds 
from a peak of X-ray emission identified by Owers et al. 
(2011). We support the general conclusion of Owers et al. 
(2011) that the major-merger in Abell 2744 is similar to that 
of the Bullet Cluster as it would be seen at a large inclina- 
tion with respect to the plane of the sky. However, we reverse 
the ordering of the major and minor mass components. 

One can infer constraints on the collisional cross-section 
of dark matter from the separation between peaks in the 
lensing and X-ray maps. For the Bullet Cluster, Markevitch 
et al. (2004) found a/m < 5 cm 2 g -1 , and for the Baby Bul- 
let, Bradac et al. (2008) found a/m < 4 cm 2 g -1 . In the Core 
of Abell 2744, we observe a projected 17" separation that, if 
the inclination is ~ 30° away from the line-of-sight (Owers 
et al. 2011), is a physical separation similar to that in the 
other bullet clusters. For an order of magnitude analysis, 
we measure the mean surface-mass density within 150 kpc 
of the mass peak £ c± 0.301q q7 g cm -2 , so that the scat- 
tering depth rs = a/m E. With the assumption r s < 1, 
which is justified due to the observed dark matter-gas sep- 
aration, we deduce a/m < 3 =b 1 cm 2 g _1 . This system may 
therefore yield one of the tightest constraints on the interac- 
tion cross-section of dark matter, based on such analysis. A 
full numerical simulation to interpret the cluster configura- 
tion would be ideal, especially given the uncertainty in the 
collision angle with respect to the plane of the sky. Indeed, 
even tighter constraints (a/m < 0.7 cm 2 g _1 ) were obtained 
from the Bullet Cluster by (Randall et al. 2008), who inter- 
preted the offsets between all three cluster components via 
tailored hydrodynamical simulations. Additional constraints 
on the collisional cross-section have also recently been ob- 



tained from dark matter stripping in Abell 3827 (Williams 
& Saha 2011; Carrasco et al. 2010) and the ellipticities of 
dark matter halos (Miralda-Escude 2002), with implications 
discussed in Feng et al. (2010). 

4.2 Northern, the bullet 

Our analysis of the Northern mass substructure (upper-left 
quadrant of Fig. 1) confirms the overall North-South merg- 
ing scenario proposed by several authors in the past. We 
find a mass ratio of ~ 2.6 between the Core and the North- 
ern clump, roughly supporting the 3:1 merging scenario of 
Boschin et al. (2006), but we identify the Northern sub- 
clump as the less massive progenitor. This conclusion is ro- 
bust, with no strong lensing features revealed by even our 
high resolution HST imaging in the Northern structure, as 
would have been expected for the reversed mass ordering 
proposed by Owers et al. (2011). 

X-ray emission in the Northern mass substructure lags 
behind the dark matter as expected. We measure a separa- 
tion of rsj 30 ;/ to the South. This is a similar separation to 
that in the core but, due to the lower surface-mass density 
in this region, constraints on the collisional cross-section are 
less significant. 



4.3 Northwestern, the ghostly and dark clumps 

By far the most interesting structure is located to the North- 
west of the cluster field (upper-right quadrant of Fig. 1). Our 
lensing analysis (Sec. 2.5) shows it to be the second-most 
massive structure. A separate region of X-ray emission also 
lies to the Northwest, called the NW interloper by Owers 
et al. (2011). Furthermore, in 54% of our weak lensing boot- 
strap realisation we identified a second peak in the more 
Western area of the NW mass clump (see Fig. 1), rendering 
it difficult to say if we indeed see a single, separate dark 
matter structure and to derive decisive separation between 
dark matter, X-ray luminous gas and bright cluster mem- 
ber galaxies. However, the separation between our different 
possible NW mass peak positions and the NW interloper is 
large (> 150 kpc for NW2). When compared to the more 
prominent Eastern mass peak NW1, which is identified in 
95% of the bootstrap realisations, the distance is even at 
least 400 kpc. 

With a de-projected temperature of ^ 5 keV (Owers 
et al. 2011) the NW interloper should have r 5 oo ~ 1-34 Mpc, 
(Evrard et al. 1996) and M 500 - 4-5 x 10 14 M o (e.g. Fig. 8 
of Khosroshahi et al. (2007)). Its 0.1 - lO.OkeV unabsorbed 
luminosity is 2.5 x 10 44 erg/s, consistent with its gas temper- 
ature of 4 — 5 keV (Khosroshahi et al. 2007; Dfaz-Gimenez 
et al. 2011). Assuming a f3 (from a King-like profile) of 0.67, 
typical for clusters, the clump would have ~ 0.95 x 1O 14 M 
within 250 kpc, similar to the N clump, and should have 
been easily detected in the lensing analysis. The interloper 
thus appears to be an X-ray feature with no associated dark 
matter or galaxies, and we therefore dub it the 'ghost' clus- 
ter. 

There is also a clear separation between the peak of 
the NW mass clump and any cluster member galaxies, so 
we call this the 'dark' cluster. Contours of the lensing mass 
reconstruction extend towards the West, where indeed we 
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find a pair of giant ellipticals (see Fig. 1). However, with the 
limited resolution of the VLT weak lensing reconstruction, it 
is impossible to tell whether this is a binary mass structure, 
though there are clear indications for that in our analysis of 
this mass clump. 

The separation between all three mass components 
makes this a real puzzle and it should be stressed that such 
a peculiar configuration is observed for the first time. It 
may pose another serious challenge for cosmological models 
of structure formation. One possible explanation was sug- 
gested by Owers et al. (2011), who describe it as a ram- 
pressure slingshot (e.g. Markevitch & Vikhlinin 2007). In 
this interpretation, after first core-passage, gas initially trails 
its associated dark matter but, while the dark matter slows 
down, the gas slingshots past it due to a combination of low 
ram-pressure stripping and adiabatic expansion and cooling, 
which enhances the cold front temperature contrast (Bialek 
et al. 2002). There is indeed a clear velocity gradient be- 
tween the NW interloper and the main cluster core (see 
Sec. 3). Such effects have also been observed e.g. in Abell 168 
(Hallman & Markevitch 2004), in a joint weak-lensing and 
X-ray study of Abell 754 (Okabe & Umetsu 2008) and in nu- 
merical simulations (Ascasibar & Markevitch 2006; Mathis 
et al. 2005), although at a much smaller separations be- 
tween dark matter and gas. The more than 100" separation 
in Abell 2744 suggests either that the slingshot scenario is 
unlikely or that some amplifying mechanism is in place. We 
shall return to this issue, proposing an interpretation of the 
entire cluster merger, at the end of this section. 

4.4 Western, the stripped clump 

The Western substructure (lower-right quadrant of Fig. 1) 
has not yet been discussed in the literature, but several clus- 
ter member galaxies are found in this area. We find a promi- 
nent weak gravitational lensing signal of ~ 1.0 x 1O 14 M0 
within 250 kpc. This should correspond to an X-ray bright- 
ness higher than the NW interloper, for the same gas tem- 
perature. However, we detect no X-ray emission. The best 
X-ray data (ACIS-I pointings 8477 & 8557) do not cover the 
region of this lensing signal, but the ACIS-S pointing 2212 
indicates no excess diffuse gas above the cosmic background 
and the extended tail of the cluster's outskirts. To match 
these observations, the Western clump must have been com- 
pletely stripped of its ICM, so we dub it the 'stripped' clump. 

The only slight excess X-ray emission nearby is a faint 
extension of the main cluster core towards the West ('ridge 
c' in Owers et al. 2011). This roughly links the Western 
clump to the Northern clump, and is consistent with rem- 
nants stripped by ram-pressure during a secondary merging 
event (NE-SW), almost perpendicular to the main merg- 
ing event projected in the sky plane. To quantify the excess 
X-ray emission, we measure counts inside three equal-size, 
non-overlapping regions extending from the cluster core to 
a radial distance of 150" ~ 680 kpc, as shown in Fig. 7. 
To prepare the data for this analysis, we divide the image 
by the exposure map following the standard procedure, re- 
move hot pixels and remove point sources using the dmf ilth 
routine in CI AO. We find a marginal excess in the central 
box (1030±32 counts), which points to the Western clump, 
above the Northern (965±31 counts) and Southern (876±29 
counts) boxes. The cluster's surface brightness profile from 
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Figure 7. Chandra ACIS-S3 X-ray image of Abell 2744, 
smoothed with a 3 pixel kernel Gaussian. We analyse the cluster's 
X-ray profile in three rectangular regions, each 150 /; in length and 
radiating from the cluster centre (marked by the blue circle). We 
find marginally significant excess X-ray emission in the central 
rectangular region, which extends towards the Western clump 
(marked by the green circle). 

the core towards the Western clump (along the central box) 
is shown in Fig. 8. The slope of the profile changes at a 
radius ~ 135" = 612 kpc away from the centre, becoming 
significantly shallower. If we take that point of transition as 
the location of the remaining gas core, we obtain a separa- 
tion between the gas and dark matter of 30". Similar results 
are found with the shallower observation 7915 using ACIS-I. 

We find a slight offset between the peak of the lens- 
ing mass reconstruction and the most luminous nearby clus- 
ter member galaxies. However, there is large uncertainty in 
the position of the lensing peak because this lies outside 
the HST imaging area. Our weak lensing analysis uses only 
VLT imaging, and there are no strong lensing constraints, so 
the mass reconstruction has a broad central plateau. Addi- 
tional HST observations would provide an ideal foundation 
to better understand the Western area, which turns out to 
be playing a significant role in the overall merger. 

4.5 One possible interpretation 

We shall now try to develop a possible explanation of the 
complex merging scenario that has taken place in Abell 2744. 
To recap, we find four mass clumps (Core, N, NW, W) with 
approximate masses 2.2, 0.8, 1.1, 1.1 x 10 14 M©, respectively. 
The Core, N and W clumps are relatively close to BCGs and 
hot gas. The NW structure, on the other hand, contains 
separated dark matter, gas and galaxies. 

We propose that the current configuration is the result 
of a near simultaneous double merger, as illustrated in Fig. 9. 
The first merger, in the NE-SW direction, had a character- 
istic path of 208" (plane of the sky distance between N and 
W clumps) or ~ 0.95 Mpc (assuming no line of sight veloc- 
ity component). The Western clump probably passed closest 
through the main cluster, as it had its ICM ram-pressure 
stripped completely. The second merger, in the SE-NW di- 
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Figure 8. Surface brightness profile (in arbitrary units) of the 
central rectangular region shown in Fig. 7. The left hand side 
starts at the X-ray centre, and the arrow denotes the approxi- 
mate position of the Western mass clump. The abrupt change in 
profile slope at roughly the same location is highlighted by the 
intersection with the horizontal line. 




Figure 9. Our proposed merging scenario, illustrated in time- 
ordered sequence, to explain the current configuration (panel 4). 
We suggest that Abell 2744 is the result of a nearly simultaneous 
double merger: one in the NE-SW direction and another in the 
NW-SE direction, which may even have consisted of three sepa- 
rate structures falling along a filament. Blue colour shall indicate 
the innermost dark matter cores of the clumps, where else their 
respective ICM is shown in red. 



rection, had a characteristic path of 117" / sin(27°) ~ 1.17 
Mpc, if we assume the inclination of that merger sug- 
gested by Owers et al. (2011) of 27°. The mergers happened 
around 0.12-0.15 Gyr ago, with a characteristic velocity of 
~ 4000 km/s as indicated by the galaxy velocity difference 
and the ICM gas velocity measurements. 

It is possible that the merger in the SE-NW direction 
could even have consisted of three initial substructures: the 
Core and two consecutive clumps (with a combined mass 
(within 250 kpc of each core) of - 1.2 x 10 14 M o ) falling 
along a filament. Those smaller clumps would be acceler- 
ated by the gravitational pull of the main cluster (plus the 
Northern and Western clumps, which were merging per- 
pendicularly). The ram-pressure slingshot in these clumps 
could be enhanced by a combination of an initially stronger 
gravitational field and perhaps a posterior reduction in ram- 
pressure due to the 'puff-up' of the gas density due to the 
recent merger of Core+N+W, similar to the density config- 
uration of the main component in the Bullet Cluster. The 
combined effect would throw the gas component ahead of its 
associated dark matter, forming the 'ghost' cluster, which is 
now the interloper. The two dark matter clumps left behind 
would now form the possibly double-centred 'dark' clump. 

This unusual scenario fits the current observations. It 
explains the clear extension of the NW clump towards the 
West and provides probably a mechanism to create such a 
huge separation between dark matter and gas in the NW 
area. However, the double peak in the NW clump needs to 
be confirmed by an additional HST pointing to cover its full 
area. The discovery of a second distinct mass peak in the 
NW clump would support our merger scenario and also de- 
crease the large measured separations between dark matter, 
gas and bright cluster galaxies. Also the ghostly interloper 
has not been covered by the current HST observations and 
the Western clump falls between chip gaps in the longer 
Chandra exposures, so further observations will still be re- 
quired. Our suggestions will also eventually require verifica- 
tion via a set of well-tailored numerical simulations and have 
to be taken with some caution at the current stage. Cosmo- 
logical boxes have to be explored (compare Meneghetti et al. 
2010, 2011) in order to find similar configurations and de- 
tailed, hydrodynamical simulations need to repeat and con- 
firm our findings. 



5 CONCLUSIONS 

We present a detailed strong lensing, weak lensing and X-ray 
analysis of the merging cluster of galaxies Abell 2744. Earlier 
studies (Kempner & David 2004; Boschin et al. 2006; Braglia 
et al. 2009; Owers et al. 2011) concluded that Abell 2744 is 
undergoing a complicated merging event. We find that it 
is even more complex than previously thought, unleashing 
a variety of exciting effects. We dub this merger therefore 
Pandora's cluster. 

Deep, three-band HST- imaging reveals a variety of 
strong-lensing features in the core of the cluster. From our 
comprehensive strong-lensing modelling of the central mass 
distribution, we identify a total number of not less than 34 
multiple images, in 11 multiple-image systems, together with 
their respective redshifts. The strong lens systems are listed 
in Tab. 1 and the finely resolved critical curve of the cluster 
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core is shown in Fig. 2. The Einstein radius of the core is 
r E = 30" ±3". 

We extended the strong lensing information with weak 
lensing measurements over the whole cluster system (600" x 
600", or no 2.7 Mpc on a side). The shapes of back- 
ground galaxies necessary for weak lensing reconstruction 
were obtained from a comprehensive combination of our new 
HST/ACS imaging, the VLT data used in Cypriano et al. 
(2004), and from additional Subaru imaging. Our combined 
strong and weak lensing mass reconstruction (Fig. 3) re- 
solves a complex structure, with at least four distinct peaks 
in the local mass distribution. The total mass of the cluster 
is 1.8 ± 0.4 x 1O 15 M within a radius of 1.3 Mpc, rendering 
Abell 2744 similar in mass to the Bullet Cluster. 

Chandra X-ray imaging also shows a complex arrange- 
ment of substructure. There are at least four X-ray ridges de- 
parting from the X-ray peak, including a Northern (NMRC) 
and a Southern (SMRC) ridge (Owers et al. 2011) . Inter- 
estingly, none of these coincide with any of the mass clumps 
found in our lensing mass reconstruction. Furthermore, the 
system also has a separate Northwestern X-ray feature with 
very low mass, also undetected in our lensing analysis. Ob- 
servations of the gas temperature in this Northwestern fea- 
ture (Kempner & David 2004; Owers et al. 2011) show a cold 
front pointing N-NE and a shock region (Tx > 15 keV) in 
the Southern ridge, for a projected impact velocity of 2150 
km/s. That projected direction (SE-NW) therefore seems to 
define the primary merging event. However, in contrast to 
the proposed merging scenario of Owers et al. (2011) , we 
find the Southern ridge to be blueshifted with respect to the 
Northern ridge. 

We also reveal that the Southern Core is ~ 2.6 times 
more massive than the Northern sub-clump - Owers et al. 
(2011) had expected this reversed - and that the secondary 
mass peak is in the Northwest. The mass ratio thus remains 
in agreement with general kinematical studies of Boschin 
et al. (2006), but swaps the sense of the collision from Owers 
et al. (2011). With this new scenario in hand, we repeated 
the velocity gradient analysis using the 'true' mass clumps 
(Core and Northwestern). These regions exhibit similar gas 
temperatures and metal abundances, and show evidence for 
a > 5200 km/s velocity gradient with the Core region red- 
shifted. These limits should be taken with caution, since the 
Chandra observations were not tailored specifically for ve- 
locity studies (they were obtained in different CCDs and at 
different epochs, so could be affected by variations in detec- 
tor gain), and there is also some evidence of unknown cali- 
bration uncertainty between the two ACIS-I pointings taken 
of the same patch of sky. However, our results are consistent 
with the main merger having a significant component along 
the line-of-sight, with a magnitude and orientation consis- 
tent with that seen in the bi-modal distribution of galaxy 
velocities. 

We also find evidence for a second merging event, simul- 
taneously with or just before the main merger. The second 
merger, along the perpendicular NE-SW axis, was between 
today's Northern and Western mass peaks. We think these 
collided inside the extended halo of the core. During this 
dramatic collision, gas in today's Northern 'bullet' clump 
was partially stripped by ram pressure, creating a charac- 
teristic separation between dark matter and baryons similar 
to that seen in the Bullet Cluster. The Western 'stripped' 



clump fared worse: all of its gas was removed, strewn into 
the tidal debris of the Core and a faint trail of excess X-ray 
emission towards its current location where we find just dark 
matter and galaxies. 

The smaller merger may have enabled a curious effect in 
the main, SE-NW merger. We postulate that gas in the Core 
was puffed up by the first collision, reducing ram-pressure 
stripping during the second. We also suggest that the main 
merger could have included two separate subclumps incident 
along a filament from the SE. The combined effect would be 
an enhancement of the 'slingshot' effect proposed by Owers 
et al. (2011), by which the subclumps' gas was accelerated 
ahead of their dark matter. This would explain the very large 
> 35" observed separation between any gas ('ghostly' clump 
or NW interloper), galaxies and dark matter ('dark' clump), 
as well as the double-peaked morphology of the 'dark' clump. 
However, this scenario needs further confirmation. 

The interpretation for this spectacular merging system 
will benefit immensely from additional observations and also 
from numerical simulations that can try to reproduce the 
new phenomenology shown in this cluster. The Western and 
Northwestern clumps have not been covered by HST obser- 
vations and Chandra observations of the Western clump are 
shallow. Wide-field Subaru imaging in good seeing condi- 
tions and in several colours would also be useful to interpret 
the global environment. Numerical simulations should be 
performed to confirm the enhancement of the baryonic sling- 
shot by the complicated merger configuration. The example 
of the Bullet Cluster has shown that the combination of 
complete lensing and X-ray observations (Markevitch et al. 
2004; Bradac et al. 2006) with highly resolved hydrodynami- 
cal simulations (Springel & Farrar 2007; Randall et al. 2008) 
is a particularly powerful tool to understand the physics of 
merging clusters. In future work, we shall attempt to repeat 
this analysis on Abell 2744. The challenge laid down will be 
to explain the complicated phenomenology associated with 
this multiple merger, as well as to better constrain the col- 
lisional cross-section of dark matter, which our rough calcu- 
lation suggests must becr/ra<3±l cm 2 g _1 , a tighter con- 
straint than that from a similar analysis of the Bullet Clus- 
ter, which does not involve numerical simulations. There- 
fore, with our own numerical follow-ups we should be able 
to place similar or even better constraints than e.g. Randall 
et al. (2008) for the Bullet Cluster. 
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